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x = 0.75, x = 0.5, and x = 0.25) are studied in a series of neutron scattering 
experiments. Powder diffraction is used to determine the temperature depen- 
dence of the magnetic structure in all four systems. Single-crystal inelastic 
neutron scattering is employed to investigate the temperature dependence of 
the Haldane-gap excitations and low-energy spin waves in the x = 1 compound 
Nd2BaNi05. The results of these experiments are discussed in the context of 
the "Haldane chain in a staggered field" model for it^BaNiOs systems, and 
quantitative agreement with theory is obtained. 



2 



Typeset using REVTgX 



I. INTRODUCTION 



The quantum-disordered ground state and the famous Haldane energy gap in the mag- 
netic excitation spectrumS have kept 1-dimensional (1-D) Heisenberg antiferromagnets (AF) 
at the center of attention of condensed matter physicists for the last 15 years. Among the 
more recent developments are studies of such systems in external magnetic fieldsiHl It was 
found that in sufficiently strong fields one of the three Haldane-gap modes undergoes a com- 
plete softening. The result is a transition to a new phase with long-range antiferromagnetic 
correlations.ifl The effect of a staggered field H^, to which a Haldane spin chain is most 
susceptible, is expected to be no less dramatic. Unfortunately, this problem has been given 
much less attention in literature, simply because such conditions are almost impossible to re- 
alize in an experiment. An arbitrarily small staggered field will induce a non-zero staggered 
magnetization in the Haldane chain. An intriguing question is how this induced order will 
affect the triplet gap excitations. Until recently the only 1-D integer-spin AF whose behavior 
was described in terms of a "staggered field" is NENP, one of the best-known Haldane-gap 
materials [Refs. |8],|9| and references therein]. Due to the existence of two non-equivalent Ni 
sites in each chain, a weak effective staggered field is induced in NENP by a uniform exter- 
nal field.0'0 While the main effect on the spin dynamics is that of the uniform field, the 
staggered component produces a small yet observable effect.0@ 

The recent discovery of coexistence of Haldane gap excitations and long-range magnetic 
order in i?2BaNi05 (i?=magnetic rare earth) compoundsEUll presents a unique opportunity 
to investigate the effect of a staggered field experimentally without having it obscured by 
the response of the system to a uniform component. According to our current understanding 
of the ii^BaNiCVs, long-range magnetic ordering does not destroy the quantum-disordered 
ground state of individual S = 1 Ni-chains in these quasi-l-D systems. It is therefore essential 
that the Ni 2+ ions carry an integer, as opposed to half-integer spin, and qualitatively differ- 
ent behavior is expected for isostructural half-integer spin materials like Nd2BaCoO5.lli0 
The non-zero staggered magnetization that appears on the intrinsically disordered S = 1 
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Ni-chains below the Neel temperature Tv is viewed as being induced by an effective stag- 
gered exchange field generated by the ordered i?-sublattice. 00 This physical picture will 
be referred to as the "Haldane Chain in a Staggered Field" (HC/SF) model for ^BaNiOs 
compounds. 

In the HC/SF model, even in the magnetically ordered state, much of the dynamic spin 
correlations are contained in Haldane gap modes propagating on the Ni chains. Below T/v 
these excitations coexist with low-energy order-parameter excitations, i.e., conventional spin 
waves, that involve correlated fluctuations of both Ni and rare earth moments. The first 
experimental confirmation of this was obtained in inelastic neutron scattering experiments 
on Pr 2 BaNiO 5 ,0 where Haldane gap modes were clearly seen both above and below the Neel 
temperature. At all temperatures these excitations have a purely 1-dimensional dynamic 
structure factor, and show no sign of softening at the 3-D AF zone-center at T = T^. 
Pr2BaNi05 samples are extremely difficult to prepare, and their magnetic properties are not 
always reproducible. For example, while magnetic ordering at T/v = 24 K has been clearly 
observed in single crystals^ for reasons yet unclear, no magnetic transition was observed 
in powders.0 It turned out that Nd-based linear-chain nickelates with the general formula 
(Nd x Y 1 _ :r )2BaNi05 are more practical as model systems for neutron studies. On one side 
of this series (x = 0) is Y2BaNiOs, a well-studied Haldane-gap S = 1 antiferromagnet 
with a spin gap A « 9 meV.&B Y 2 BaNi05 remains magnetically disordered even at low 
temperatures. At the other end we have Nd2BaNiOs (x = 1) that, thanks to the presence of 
magnetic Nd 3+ ions, orders antiferromagnetically at T/v = 48 K. As in the case of Pr 2 BaNi05, 
Haldane-gap modes in Nd 2 BaNiOs were found to survive in the magnetically ordered phase.0 
Similar behaviour was seen in x = 0.75, x = 0.5 and x = 0.25 species, where the ordering 
temperatures are 39 K, 30 K, and 19 K, respectively^ For all Nd concentrations at T > T N 
the Haldane gap modes are virtually indistinguishable from those seen in Y 2 BaNiOs. As the 
temperature is decreased through T/v however, the gap energy A increases roughly linearly 
with T/v — T. This remarkable behaviour is in agreement with theoretical predictions based 
on HC/SF modelS! 
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To date, all inelastic measurements on (NdaYi-^BaNiOs materials were performed on 
powder samples. The existing magnetic structure (diffraction) data are not accurate or not 
complete enough to be discussed in relation to our model. Indeed, previous single-crystal 
experiments on the Nd-system@0 were analyzed assuming that the Nd and Ni moments 
are collinear at all temperatures below T^. It was later shown that the actual structure 
is canted,il but no detailed temperature dependencies derived from the canted model were 
reported. For the mixed Nd-Y systems only very preliminary results obtained by analyzing 
only three magnetic reflections are available .0 

In this work we continue our studies of (NdaYi-^BaNiOs compounds using both elas- 
tic and inelastic neutron scattering. Powder neutron diffraction is used to investigate the 
temperature evolution of the magnetic structure in four (Nd^Yi^^BaNiOs species with 
x = 1, x = 0.75, x = 0.5 and x = 0.25, respectively. Inelastic neutron scattering is then 
employed to study the temperature dependence of the Haldane-gap excitations in the Ni- 
chains as well as the low-energy order-parameter fluctuations (spin waves) in single-crystal 
Nd2BaNiC>5 samples. The connection between static and dynamic properties is discussed 
in the framework of the HC/SF model. Some of the results presented below were briefly 
discussed in Ref. I2T5. 



II. EXPERIMENTAL PROCEDURES 

The crystal structure of i? 2 BaNi0 5 compounds is discussed in great detail elsewhere.0 
All (Nd x Y 1 _ 2 .) 2 BaNi05 systems crystallize in an orthorhombic unit cell, space group Immm, 
a ~ 3.8 A, b ~ 5.8 A, and c ~ 1 1. 7 A.ll! The S = 1 Ni chains are formed by orthorhombically 
distorted Ni06 octahedra that are lined up in the [100] crystallographic direction, sharing 
their apical oxygen atoms. The R sites (occupied at random by Y or Nd) are positioned in 
between these chains. The site symmetry for the R 3+ is low and the electronic ground state 
for Nd 3+ is a Kramers doublet. 

In powder diffraction experiments on the x = 0.75, x = 0.5 and x = 0.25 compounds 
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we used the same samples (roughly 15 g each) as in previous inelastic studies!! A new 
Nd2BaNiC>5 (x — 1) powder sample of comparable mass was prepared particularly for the 
present study using the solid state reaction method. Powder neutron diffraction experiments 
were carried out at the high flux reactor of the Institut Laue-Langevin (Grenoble), on the 
400 cells position sensitive detector diffractometer DIB for the x = 0.75, x = 0.5 and 
x = 0.25 compounds and on the 1600 cells diffractometer D20 for the x = 1 system. The 
patterns were recorded in the temperature range 2 - 50 K, using a liquid helium cryostat. 
The wavelengths A = 2.526 A on DIB and A = 2.400 A on D20 were provided in both cases 
by focusing pyrolitic graphite monochromators. The samples were enclosed in cylindrical 
vanadium containers. The analysis of the powder patterns was performed by Rietveld profile 
refinement using the software FULLPROF.iil A pseudo-Voigt function was chosen to generate 



the line shape of the diffraction peaks. The scattering lengths were taken from Ref. 32 and 



the magnetic form factors of Nd and Ni from Ref. |33 

Single-crystal samples for inelastic neutron scattering experiments were grown in oxy- 
gen atmosphere using the floating-zone technique. Most of the thermal neutron inelastic 
measurements were done on a cylindrical single crystal 25 mm long and 4 mm in diameter 
(sample A) with mosaic about 0.3°. Only near the end of this series of experiments did a 
much larger sample (sample B) become available. Sample B consists of three single crystals, 
each roughly 6 mm in diameter and 50 mm long with an « 0.8° mosaic. The three crys- 
tals were aligned together on a specially designed aluminum sample holder. The resulting 
"supersample" had a symmetric mosaic of around 1.2°. Sample B was used primarily in 
cold-neutron inelastic scattering experiments. 

Inelastic thermal neutron scattering measurements were carried out at the Center for 
Neutron Research at the National Institute of Standards and Technology (NIST) on the 
BT-4 and BT-2 triple-axis spectrometers. The samples were in all cases mounted with the 
(h, 0, 1) reciprocal-space plane parallel to the scattering plane of the spectrometer. A neutron 
beam of 14.7 meV fixed final energy was used with 60' — 40' — 40' — 120' collimations and a 
pyrolitic graphite (PG) filter positioned after the sample. PG (002) reflections were used for 



monochromator and analyzer (setup I). The measurements were done with energy transfers 
up to 25 meV, the calculated resolution at a typical 15 meV energy transfer being 2.5 meV 
FWHM. The sample environment was a standard displex refrigerator, and measurements 
were done in the temperature range 25-60 K. As we have found, in all samples the axis 
of the cylindrical crystal roughly coincides with [100] crystallographic direction. In the 
experiment the sample cylinder was therefore positioned horizontally in the scattering plane, 
which resulted in substantial absorbtion effects. The effective transmission coefficient as a 
function of sample orientation was determined by measuring the intensity of incoherent 
scattering. The measured absorption correction was almost constant within each constant- 
Q scan performed, but varied by as much as a factor of two between different Brillouin 
zones. 

Inelastic cold neutron scattering experiments were done on the SPINS 3-axis spectrome- 
ter at NIST. Several configurations were used. In setup II we utilized 30' — 80' — 80' — open 
collimations with a 3.7 meV fixed final energy neutron beam and a BeO filter after the sam- 
ple. In this setup the energy resolution at zero energy transfer is 0.2 meV FWHM. 5 meV 
fixed-final energy neutrons were used in setup III with the same collimation setup and a 
Be-filter positioned after the sample. A substantial gain in intensity (at the cost of wave 
vector resolution) was obtained using 30' — 80' — 45'(radial) — open collimations and 5.1 meV 
fixed final energy neutrons, a Be-filter after the sample and a horizontally-focusing analyzer 
(setup IV). The energy resolution in this case was 0.35 meV FWHM at zero energy transfer. 
The sample environment was a displex refrigerator and the temperature was controlled in 
the range 10-50 K. 



III. EXPERIMENTAL RESULTS 
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A. Magnetic structure 



For all samples studied the diffraction patterns in the paramagnetic phase were found to 



be totally consistent with the crystal structure reported in Ref. 30. A small amount of non 
magnetic impurities was detected in the x = 0.75, x = 0.5 and x = 0.25 samples, whereas 
the x = 1 one turned out to be single phase. The small extra peaks due to this impurity 
were excluded from the refinements. 

As previously observed^ magnetic ordering in all samples manifests itself in the ap- 
pearance of new Bragg reflections at half-integer positions, the magnetic propagation vector 
being [1/2 1/2]. In the x = 1 sample, thanks to the relatively strong magnetic signal, a 
simultaneous refinement of nuclear and magnetic structures could be performed. A typical 
powder scan for this system is shown in Fig.|I](a). For the other three compounds, where 
the magnetic signal was weaker, a different approach was used. For powder scans collected 
at T < Tjv the magnetic contribution to scattering was isolated by subtracting the nu- 
clear background measured just above the ordering temperature. Several typical data sets 
obtained in this fashion are shown in Fig.[|(b-d). 

In the x = 0.25 sample, in addition to the above mentioned magnetic peaks at half- 
integer positions, extra intensity below T/v was also observed at the positions of nuclear 
peaks [Figj^(d)]. One possible explanation for this is the presence of an extra magnetic 
component in the structure of this compound. This however appears not to be the case. We 
were unable to reproduce the observed pattern with any reasonable spin model. In addition, 
the intensities of these extra reflections seem to increase with increasing scattering angle, 
strongly suggesting that they are of nuclear, rather than magnetic origin. In our analysis 
we have assumed that theses extra peaks represent a lattice distortion that is induced by 
magnetic ordering. In the determination of the magnetic structure of this sample these 
reflections were therefore ignored. 

According to the analysis of the crystallographic space group of R^BaNiOs compounds 
(Immm), only two types of magnetic structure with propagation vector [1/2 1/2] and 
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non-zero magnetic moment on both Ni and R sites are possible.Ej In the structure realized 
in Er 2 BaNi05 all spins are confined to the (a, b) crystallographic plane. As shown by several 
previous studies,@@ it is the other alternative that is realized in Nd2BaNiOs: both the Ni 2+ 
and Nd 3+ moments are in the (a, c) plane of the crystal (Fig. |2|). We used this model to fit the 
diffraction spectra measured at each temperature for each sample. The adjustable parameters 
were the magnetic moments m*^ 1 ' and m^ Nd - ) of the Ni 2+ and Nd 3+ ions, respectively, as well 
as the angles <^ Nl ) and <^ Nd ) these moments form with the c axis of the crystal. Good fits 
to the experimental powder profiles, examples of which are shown in solid lines in Fig. [l], 
were obtained in all cases. The refined values of parameters are plotted against temperature 
in Fig. |3] and Fig. f|. As observed in previous studies,^ near the Neel temperature the 
magnetizations of both Ni and Nd sublattices undergo a slight reorientation with decreasing 
temperature. Very quickly though this reorientation slows down and in all samples the angles 
0( Nl ) anc [ 0( Nd ) remain at roughly —35° and 0°, respectively, in most of the temperature range. 



B. Ni-chain gap excitations 

The temperature dependencies of energies and intensities of the Haldane-gap excita- 
tions were studied by inelastic thermal neutron scattering. One of our objectives was to 
use the intrinsic polarization dependence of the neutron scattering cross section to distin- 
guish between the three individual modes, expected to be present in the Haldane multiplet. 
For this purpose the data were collected in constant-Q scans at the 1-D antiferromagnetic 
zone-centers Q = (1.5,0,0) and Q = (0.5,0,4.2). These wave vectors correspond to equal 
momentum transfers, but are pointing along and almost perpendicular to the Ni-chain axis 
a, respectively. At (1.5, 0, 0) only the fluctuations of y and z spin components are seen (the 
axes x, y and z are chosen along the [100], [010] and [001] directions, respectively), while at 
(0.5,0,4.2) it is mostly the fluctuations of x and y spin components that contribute to the 
scattering intensity. 

The most serious obstacle in previous powder measurements was the presence of two crys- 
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tal field (CF) excitations associated with Nd 3+ that appear at about 18 meV and 24 meV, 
respectively.lll Unlike the highly dispersive Haldane excitations, the dispersionless CF modes 
do not suffer an intensity penalty upon powder averaging the corresponding dynamic struc- 
ture factors. The Haldane modes that appear at « 10 meV energy transfer at T = 50 K 
and move towards higher energies at lower temperature are thus very difficult to separate 
from the intense CF background. In a single crystal this problem may be overcome. First, 
the intensity of Haldane modes relative to those of CF excitations is much larger in this 
case. Second, in a single crystal the background can be directly measured at scattering 
vectors slightly off the 1-D zone-center, where the Haldane excitations move out of the 
scan range thanks to their steep dispersion. In our case the background was measured at 
Q = (1.4,0,0) and Q = (0.4,0,4.2). Since it was previously found that the CF modes are 
slightly temperature-dependent, the background measurements were performed separately 
at each temperature. Figure [5] shows some typical data collected in sample A at T = 50 K 
(open circles) together with the corresponding background scans (solid circles). The result 
of point-by-point subtraction is shown in Fig. |^. The well-defined peak at ~ 11 meV en- 
ergy transfer, very similar to that previously seen in Y 2 BaNi0 5 iHI and Pr 2 BaNi05,t3 is 
attributed to the Haldane excitations in the Ni chains. 

One problem that had to be dealt with was the presence of a spurious peak at roughly 
18 meV energy transfer (see Fig. |^). We have identified this spurion as being caused by 
higher-order scattering (2&j) in the monochromator, incoherent and/or Bragg scattering in 
the sample, and higher-order (3k f) scattering in the analyzer. In the analysis described 
below the "dangerous" energy range has been excluded from all data sets. 

Ideally, for a complete determination of the mode polarization, at each temperature 
one would analyze the inelastic scans collected at both wave vectors simultaneously, using a 
model cross section that would include three modes with separate gap values and polarization 
factors. In the present study however, as the structure of the observed peak is totally smeared 
by the broad instrumental resolution, using a three-component model cross section results 
in an over-parameterized problem. Instead we analyzed all the measured inelastic scans 
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separately using a single- mode cross-section, as was previously done for P^BaNiOs. B The 
model dynamic structure factor was written in the "double-Lorentzian" form: 

S(q,u) = - 77— V2TF2 i 1 ) 

(^ q ) 2 = A 2 + c 2 q 2 . (2) 

In this formula q is the projection of the scattering vector q onto the chain axis, and measured 
relative to the 1-D antiferromagnetic zone-center; A is the gap energy; c is the spin wave 
velocity along the chains; T is the intrinsic energy width of the excitations; and £ is the real- 
space dynamic spin correlation length. Since all measurements were performed at q — 0, 
the parameters £ and c are only needed to properly take into account focusing effects. In 
our analysis £ and c were therefore fixed to the values previously measured in P^BaNiOs: 
1 /£ = 0.08 A -1 , and c = 200 meV A. For each scan the cross section ([!]) was folded with the 
4-dimensional spectrometer resolution function and the parameters A, V and the prefactor 
So were refined to best-fit the data. To further minimize the number of independent variables 
we have averaged the value of T determined at different temperatures and used this average 
as a fixed parameter in the final fit. This was done separately for Q = (1.5, 0, 0) and 
Q = (0.5,0,4.2), where the average values were 2 meV and 1.5 meV, respectively. Several 
typical fits of Eqs. to our inelastic data are shown in solid lines Fig. 0. The obtained 
temperature dependencies for A and So for both wave vectors are shown in Fig. [|. 



C. Low-energy spin waves 

The main purpose of the cold-neutron experiments was to study the low-energy order- 
parameter fluctuations in Nd2BaNiOs. Usually, such excitations resulting from long- 
range ordering appear as acoustic spin waves, as previously observed in Pr 2 BaNiO 5 .0 For 
Nd2BaNiC>5 most of the measurements were done in the vicinity of the (0.5, 0, 1.5) magnetic 
Bragg reflection, one of the strongest magnetic peaks in the (a, c) plane, with roughly 10000 
counts per second at T = 10 K using setup II. To search for the acoustic spin waves we 
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used configurations I, II or III to perform constant-Q scans at wave vectors (0.5,0,1.5), 
(0.75, 0, 1.5), (1, 0, 1.5), and (1.5, 0, 1.25) in the range 0-8 meV (Fig. g), as well as constant- 
E scans at an energy transfer Too = 0.3 meV along (0.5 + £,0,1.5) and (0.5,0,1.5 + £) 
(—0.25 < C < 0.25). At T = 10 K absolutely no inelastic signal was found at energy 
transfers below rs 4 meV. At hu = 4 meV we observe the rather intense resolution-limited 
inelastic feature previously seen in powder experiments.0 The 4 meV excitation is totally 
dispersionless along both the a and c axes of the crystal. This can be deduced from Fig. [H| 
that shows some typical constant-Q scans collected at several wave vectors at T = 10 K. 
Scanning through the 4 meV peak in different Brillouin zones led us to the conclusion 
that the mode is a single-ion excitation with no unit-cell structure-factor and no apparent 
polarization-dependence of intensity. The accessible range of momentum transfers was not 
sufficient for an accurate study of the form-factor dependence of the intensity. 

What makes the 4 meV mode rather interesting is its temperature dependence. As was 
previously established in inelastic nentron scattering experiments on powder samples.B this 
excitation is visible only at T < Tjv- As the Neel temperature is approached from below the 
excitation energy decreases, the width increases dramatically (Fig.[TT|), and at T = T N the 
mode appears to merge with quasielatic scattering at Tiu = 0. We have utilized setups I and 
IV to measure the temperature dependence of the excitation energy for temperatures up to 
T = 35 K, where the inelastic peak is broad, yet still easily identified. Constant-Q scans 
were analyzed with a damped-oscillator dynamic structure factor: 

1 1 



_ Iq7 [n(uj) + 1] 
b{u — 

UJqTT 



(3) 



that was convoluted with the gaussian energy-resolution of the spectrometer (solid lines in 
Fig.|TT|). In Eq. ^| n{oj) is the Bose factor and I is the structure factor of the excitation. 
For Q = (0.5,0,1.5) the experimental T-dependencies of the excitation energy Tiujq and 
relaxation rate 7 are shown in Fig. [12|. The same behavior was also seen at Q = (0.75, 0, 1.5) 
and Q = (1,0,1.5). 
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IV. DISCUSSION 



A. Magnetic ordering 

In previous resonant magnetic X-ray diffraction studies^ it was found that the order- 
parameter critical exponent (3 in N^BaNiOs is indistinguishable from the mean-field (MF) 
value /3 = 0.5. In fact, for the x=l system, the standard MF theory can accurately describe 
the entire temperature dependence of both Ni and Nd ordered moments.0 In this model 
both Ni 2+ and Nd 3+ ions are treated as isolated spins in an effective exchange field, and 
their bare magnetization curves are given by appropriate Brillouin functions. The reduced 
moment on the Ni-sites (the expected value is 2 /x B ) is attributed to orbital effects (reduced 
g-factor). Despite that reasonably good fits to the experimental temperature dependencies 
can be obtained by applying this simple approach to each N^BaNiOs compound separately, 
the refined Ni-Nd and Ni-Ni coupling parameters, as well as the ^-factor for Ni 2+ depend 
greatly upon sample composition. Such a strong variation of the Ni-Ni exchange constant is 
difficult to explain by the slight composition dependence of the lattice parameters and hence 
the Ni-Ni exchange integral.0 Indeed, the Haldane gap energies in the paramagnetic phase 
are directly proportional to the Ni-Ni exchange constant, and were previously found to be 
independent of Nd content .El If we are to describe the measured temperature dependencies 
of sublattice magnetizations in (Nd^Yi-^BaNiOs, treating magnetic interactions at the 
MF level, we have to look for a different approach. 

In-chain Ni-Ni AF exchange coupling is by far the strongest magnetic interaction in 
i? 2 BaNi05 compounds (Jm-ni ~ 300 K). At T < A, which, incidentally is the temperature 
range where long-range magnetic ordering occurs, the susceptibility of the Ni subsystem 
is therefore dominated by well-established in-chain correlations. In this situation it is not 
appropriate to model the magnetic transition by treating the Ni 2+ ions as isolated moments. 
Instead, the MF theory should treat entire Ni-chains as inseparable entities, characterized 
by their bare staggered magnetization function. 3-dimensional long-range ordering is driven 
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by the next important magnetic interaction, namely Ni-R exchange coupling that establishes 
links between individual chains and is of the order of tens of Kelvin. Direct R-R interactions 
are expected to be very weak (of the order of 1 K or less) and can for this reason be 
disregarded.^ 

Garcia-Matres et al.0 were probably the first who analyzed their magnetic susceptibility 
data treating only the Ni-R interactions at the MF level. The experimentally measured 
temperature dependence of the Ni-order parameter was built into this model as an empirical 
function. The HC/SF model provides a physical justification for this approach:!! the bare 
magnetization curve for the Ni-subsystem is to be taken in the form of the staggered magne- 
tization function for an isolated quantum-disordered Haldane chain. This "semi-quantum" 
model differs qualitatively from the standard MF model for classical magnets. In the latter, 
at T — > all sublattices become fully saturated as the bare single-ion susceptibility diverges 
as 1/T. In contrast, the bare staggered susceptibility Xn of a Haldane spin chain remains 
finite at T = 0, and the Ni-sublattice need not be fully saturated. The "semi-quantum" 
model can thus account for the reduced ordered moment on the Ni-sites without assuming 
an improbably small g-factor for Ni 2+ . Moreover, it will predict qualitatively different tem- 
perature dependencies, since the staggered susceptibility of a Haldane spin chain, unlike the 
single-ion susceptibility, is expected to be almost T-independent at T < A (A ps 9 meV, or 
pa 105 K in i? 2 BaNiOs compounds). 

It has to be noted that treating interchain coupling at the MF or RPA (Random Phase 
Approximation) level is, in itself, not a new idea, but a well-established technique!!. In 
their pioneering work on CsNiCls, the first Haldane-gap material studied experimentally, 
Buyers et al.@, and Affleckill implemented this approach to explain magnetic ordering and 
calculate the spin wave dispersion relations. The main difference between CsNiCls and 
i? 2 BaNi05 is that exchange coupling between individual Haldane chains is direct in the 
former system, and mediated by the rare earth ions in the latter. For directly coupled 
Haldane spin chains, the magnitude of interchain interactions must exceed some critical 
value in order for the system to order magnetically^. In a MF description of CsNiCls or 
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any other system with directly couple quantum chains, the temperature dependence of the 
ordered moment is defined by the intrinsic temperature dependence of the susceptibility of 
individual chains. In contrast, in our case of ^BaNiOs compounds, magnetic ordering is 
driven by the 1/T-divergent susceptibility of the rare earth subsystem. At sufficiently low 
temperature long-range order will therefore occur for arbitrary small R-Ni interactions, and 
in the case of T/v A, we can use the approximation where the bare magnetization curve 
of isolated chains is T- independent. 

An important results of this work is that in all samples studied , the induced stag- 
gered cmoment on the Ni-chains, appears to be a universal function of m^ Nd ^ = m^^x. 
In other words, the Ni magnetization explicitly depends only on the average moment on 
the i?-sites, and not the actual temperature, as illustrated in the plot of m^ Nl ) vs. m( Nd ) 
in Fig. [13|. To emphasize the significance of this fact we reiterate that the ordering tem- 
peratures in the materials differ by more than a factor of two. From Fig. we see that 
even with fully ordered Nd 3+ moments the Ni-chains are not fully polarized and the curve 
never levels off completely. In the context of the HC/SF model, apart from the scaling of 
the absciss, Fig. [TB| is nothing else but the staggered magnetization function Ai(H n ) for an 
isolated Haldane chain of Ni 2+ spins, that in the studied temperature range T = — 50 K 
is expected to be almost T- independent. 

The role of disorder in Y-diluted systems is not too important, at least, in systems where 
the Nd concentration is not too small. The intrinsic dynamic spin correlation length of 
Haldane spin chains is rather large, of the order of six lattice repeats. As long as the mean 
distance between the Nd 3+ ions is smaller than this length scale, the Ni-chains effectively 
see a homogeneous staggered moment on the i?-sites, despite the "holes" that are present 
wherever Nd 3+ is replaced by Y 3+ . The argument, while clearly valid at T « T/v, should be 
taken with some caution when applied to a system deep in the ordered phase. In this regime 
the Haldane gap energy increases, and, as a consequence of that the dynamic correlation 
length in the chains decreases. Even in this case we can expect the effect of disorder to be 
averaged out, thanks to the simple geometrical fact that every Ni 2+ ion is coordinated to 
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four .R-sites. 

We can now use the experimentally determined staggered magnetization function for 
Haldane spin chains in (Nd^Yi-z^BaNiOs to write down the self-consistent MF equations 
for our "semi-quantum" theory. For the Ni-sublattice these equations are: 



m 



(Ni) 



gSfx B M(H^) (4) 
= 2cW Nd > (5) 

Here and S — 1 are the gyromagnetic ratio and spin of Ni 2+ , respectively; if( Nl ) is the 
effective exchange field that acts on the Ni-chains and is generated by the i?-sublattice; and 
a is an effective MF coupling constant. The factor 2 in Eq. |5| reflects the fact that there are 
two Nd atoms for every Ni atom in the chemical formula, while the magnetizations mP* 1 ' 
and m ( - Nd ^ are normalized per site. For convenience .M(-£^ Nl )) in Eq. |] is approximated with 
the following purely empirical function: 

gSfi B M{2am {Nd) ) = A arctan(Em (Nd) ). (6) 

The coefficients A = 1.27(5) and B = 1.06(13) are obtained by fitting this formula to the 
bulk of experimental data in Fig. [H| (solid line). To write down the remaining MF equations 
for the i?-sublattice we shall approximate the bare magnetization curve for the Nd 3+ ions 
with the Brillouin function, as is appropriate for a magnetic ion with a doublet ground state: 

m (Nd) = m (Nd) tanh( 4 Nd) ^ (Nd) ) (?) 

kT 

H^ = am^. (8) 

In this formula if ( Nd ) is the effective exchange field acting on the rare earth ions, n is 
Boltzmann's constant and m!^^ is the effective magnetic moment for Nd 3+ . 

Ideally, the saturation moment mo should be indepent of temperature or Nd concentra- 
tion in the sample. In our particular systems however, the saturation moment of the Nd 3+ 
ions steadily decreases by roughly 25% as x changes from 1.0 to 0.25. Eq. |7| is therefore no 
more than a crude approximation, and a 25% accuracy is the best we can expect from our 
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model. Indeed, in using the Brillouin function to describe the rare earths, we have totally 
neglected the higher-energy electronic states of these ions. Not only do they contribute to 
the temperature dependence of i?-magnetization, but also may give rise to a T-independent 
Van-Vleck contribution to single-ion susceptibility. Unfortunately, without knowing the elec- 
tronic structure of R 3+ in R 2 BaNi0 5 in detail, we can not take these effects into account 
rigorously. Instead we shall use Eq. [7| with separate values for mo for each sample, equal to 
the saturation magnetization of the Nd sublattice. 

Equations ([§-§) were used to analyze the temperature dependencies of sublattice mag- 
netizations in three (Nd^Yi-^BaNiOs samples, with x — 1, x — 0.75 and x = 0.5. In the 
x = 0.25 system the transition is of magneto-elastic nature, and a purely magnetic MF model 
can not be expected to be applicable. For each sample composition the only adjustable pa- 
rameter was the coupling constant a. The results of this analysis are shown in solid lines in 
Fig. g. For the three samples we obtain a = 0.84(2) • 10 5 Oe /i^ 1 , a = 0.96(2) • 10 5 Oe fig 1 , 
and a = 1.14(2) ■ 10 5 Oe , respectively. As expected, for all three (Nd^Yi^^BaNiOs 
systems the refined values of the coupling constant are very similar despite the substantially 
different ordering temperatures. We can use average value a = 0.98(1) • 10 5 Oe /i^ 1 to 
properly rescale the absciss in Fig. [T3] and obtain the Ni-order parameter as a function of ef- 
fective exchange field, i.e., the actual staggered magnetization curve for a Haldane spin chain 
(Fig. [13], top axis). Differentiating Eq. ^| at — > we obtain the staggered susceptibility 
X, = 1.2 • 10- 5 /VOe. 

The staggered magnetization curve for an isolated S = 1 Heisenberg AF chain has not 
been calculated to date, and we can not directly compare our results to any numerical 
simulations. The zero-temperature staggered susceptibility Xn, on the other hand, can 
be deduced from the numerous theoretical predictions for the dynamic structure factor 
at the 1-D AF zone-center S(tc,uj), to which it is related through the Kramers-Kronig 
relation and the fluctuation-dissipation theorem. In the single-mode approximationilSi 
S(ir,u) = — ^^8(uj — A). Here and ^ is the ground state energy per spin. According 
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to Monte Carlo simulations (see, for example, Refs. gggl]), ^ « -1.4J « 3.42A. The 
expression for the staggered susceptibility is then readily obtained via x-k = 2 J °° du S ^' w ^ 
and is equal to: 

Xw = 9.1(^ fl ) 2 /A. (9) 

In all the i? 2 BaNi05 systems studied so far with inelastic neutron scattering the gap energies 
are equal to 9 meV to within experimental error, which gives x-k ~ 2.3 • 10 _5 /is/Oe. This 
value is almost twice as large as our experimental estimate. While not spectacular, this level 
of consistency is quite acceptable, considering all the approximations and simplifications that 
had to be made. Most of the uncertainty is related to the large experimental error in the 
initial "low-field" part of the magnetization curve in Fig. 12 where small ordered moments 
were derived from the measurements in the vicinity of T/v, and the error bars are rather large. 
Some systematic error is introduced by approximating the measured curve with the expres- 
sion (H). We have also totally neglected the temperature dependence of the magnetization 
function for the Haldane spin chains, small as it may be at kT < 0.5A, and approximated 
that for the rare earth ions by a simple Brillouin function. The actual gyromagnetic ratio 
for Ni 2+ that enters Eq. || squared, is not known from independent measurements either. A 
half-order-of- magnitude agreement is indeed the most we can expect. What is important 
is that our model appears to be self-consistent and is based on the same concepts as those 
previously used to explain the persistence and T-dependence of Haldane-gap excitations in 
i? 2 BaNi05 materials. 



B. The 4 meV mode 

The excitation that we see at 4 meV energy transfer was first observed in powder 
experiments,0 and attributed to a CF transition in Nd 3+ . A similar feature was also seen in 
recent AFMR measurements.^ As our new single-crystal neutron data clearly demonstrate, 
the 4 meV mode has indeed no dispersion and its intensity is Brillouin-zone independent, 
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so it indeed is a local, i.e., single-ion excitation(s). Despite that, we suggest that this mode 
is the actual order-parameter excitation (spin wave), and corresponds to flipping a single 
Nd 3+ moment in the effective staggered field projected by the Ni-sublattice. The energy of 
such an excitation is given by: 

hu = 2m # (Nd) = 2m am (Ni) . (10) 

The excitation energy should thus scale as the ordered moment of the Ni-sublattice. The 
solid line in Fig. |12| is drawn using Eq. [U|, our experimental data for m^ Nl -*(T) in Nd2BaNiOs, 
and a = 0.84 ■ 10 5 Oe/x^ 1 determined from our MF analysis of the magnetic structure. An 
excellent agreement with inelastic measurements is apparent. 

One possible explanation for the fully localized nature of the low-energy spin waves is a 
strong anisotropy in the effective exchange interaction between Ni 2+ and Nd 3+ . We suggest 
that this interaction is of Ising type, i.e., can be written as Je«SmS§ideSi SNd.efr being the 
effective Nd-spin. If one introduces this type of coupling between singlet ground state Ni- 
chains and isolated paramagnetic rare earth ions at the RPA level, one obtains a total 
decoupling of the Ni- and Nd- transverse degrees of freedom. The Ni-chain modes are those 
of a Haldane system in a static staggered field. Those centered on Nd correspond to flipping 
a single paramagnetic moment, again in a static effective field. Interestingly, the low-energy 
excitations in Pr 2 BaNi05 are qualitatively different: these are acoustic spin waves with a 
pronounced dispersion and structure-factor .111 These excitations, as well as the higher-energy 
Haldane-gap modes, are expected to involve correlated fluctuations of both the Ni and the 
rare earth moments. The difference in behavior is, of course, due to the different electronic 
configuration of the rare earth ions in P^BaNiOs and Nd2BaNiOs (Pr is not a Kramers ion 
and in Pr 2 BaNiOs is, in fact, an induced-moment system), and the consequent difference in 
R-Ni magnetic interactions. 

Unfortunately, within the MF model one can not make any predictions regarding the 
strongly T-dependent damping of the 4 meV mode. Without going into further speculations, 
we would only like to point out that the measured magnitude and temperature dependence of 
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7 is very similar to that of the temperature-induced broadening of Haldane-gap excitations, 
as measured in Y 2 BaNi0 5 .H One possibility is that these quantities may actually be directly 
related and represent the same relaxation process. 



C. Haldane gap excitations in the Ni chains 

1. Polarization 

One set of results that we hoped to obtain from the inelastic thermal neutron scattering 
experiments are separate temperature dependencies of the Haldane gap energies in each of 
the three components of the Haldane triplet. Of particular interest is the T-dependence of 
the gap in the longitudinal mode, i.e., the one polarized along the ordered moments on the 
Ni-sites. It has been predicted that the rate of increase of the longitudinal spin gap is three 
times larger than that of the two transverse gaps.@. As mentioned above, experimental 
difficulties make a complete polarization analysis impossible. At present, the only available 
results relevant to the multiplicity and polarization of Ni-chain gap excitations in Nd2BaNiOs 
is the temperature dependence of integrated intensities measured at Q = (1.5,0,0) and 
Q = (0.5,0,4.2) (Fig. [5J). To within experimental error these intensities are the same at 
all temperatures. At no instance do we see any distinct splitting in the energy gap. If we 
assume that the modes remain polarized along the principal crystallographic directions, as 
dictated by the local symmetry of the Ni-sites, we are forced to conclude that the three 
components of the triplet evolve with temperature in exactly the same way. 

The assumption however is expected to be valid only at T > T^. As the magnetic 
order-induced shift in the Haldane gap energies becomes larger than the initial anisotropy- 



induced splitting of the triplet (roughly 1 meV, according to Refs. |24l , |2~3|) , the proper choice 
of polarization axes is governed by the direction of the Ni 2+ moments, rather than the 
crystallographic symmetry directions. In the magnetically ordered state the three excitations 
are polarized along the x', y', and z' axes, where y' runs along the [010] direction, z' is chosen 



20 



along the direction of ordered moments on the Ni sites and thus forms an angle of roughly 
35° with the c-axis of the crystal, and x' complete the orthogonal reper, respectively. The 
actual situation is even more complex, since two types of magnetic domains with the Ni 2+ 
moments tilted to the left or to the right of the c-axis should be present in a macroscopic 
sample. As a result all three modes contribute to inelastic scattering at both Q = (1.5, 0, 0) 
and Q = (0.5, 0, 4.2), with partial intensities 0.33, 1, 0.67, and 0.63, 1, 0.37, correspondingly. 
We have concluded that our data are not inconsistent with seeing only the two transverse 
excitations (polarized along x' and z', respectively). It is therefore entirely possible that for 
some reason we fail to see the longitudinal excitation alltogether. 

Definitive conclusions concerning the multiplicity and polarization of the gap excitations 
are premature. At this stage however, we lean towards interpreting our data in terms of 
seeing only the transverse Haldane excitations in our scans. First of all, the longitudinal 
gap is expected to increase with decreasing T at least three times as fact as the transverse 
gaps. Below T = 40 K it is expected to be already outside the range of our inelastic scans. 
Second, the observed temperature dependence of the Haldane gap energies is in striking 
quantitative agreement with theoretical predictions for transverse modes (see next section). 
Third, it appears plausible that the appearance of an ordered moment on the Ni-sites leads 
to the longitudinal mode being overdamped. Hopefully, the use of the very large sample B 
in future unpolarized and polarized neutron scattering experiments will help us resolve this 
issue. 

D. Behavior of the Haldane gap energy 

As discussed in the previous section, our results for the low-energy spin waves indicate 
that magnetic excitations involving the Ni-spins are totally decoupled from those in the 
Nd-subsytem and are indistinguishable from transverse Haldane gap excitations that occur 
in isolated Haldane chains in the presence of a static staggered field. The HC/SF model can 
thus be used in reference to spin dynamics, and not only the static magnetic properties of the 
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system. Based on this assumption (that, at the time, was not backed by experimental results 
pertaining to the low-energy spin waves in the system), the increase of the Haldane gap 
energy in the magnetically ordered state in (NdzYi.^BaNiOs was qualitatively explained 



in Ref. [18] In this work the staggered-field dependence of the Haldane gap energies for a S = 1 
1-D Heisenber antiferromagnet was derived through the use the use of a Ginsburg-Landau- 
like description of an isolated Haldane spin chain, the so-called 4 -model, first introduced 
in this context by Ameck.0il 

In this approach an isolated Haldane spin chain is characterized by the space-dependent 
vector of local staggered magnetization <p(x,t). The Lagrangian of the system is written as 



a series expansion in <p{x) and its derivatives: 



£ = dx 



^("Ir) 2 ~ Stir") 2 ~ TT<f> 2 ~ + higher-order terms 
2v at 2 ox 2v 



(11) 



The quadratic terms describe a triplet of degenerate non-interacting gap excitations (energy 
gap A) with a linear dispersion (spin wave velocity v ) . A non-zero staggered-field dependence 
of the excitation energies appears upon introducing the A|0| 4 -term, that represents pair- 
repulsion between magnons, and other higher-order terms. In the presence of a staggered 
field, the change in gap energies, defined as S(H W ) = ^ H *A — A being the gap energy in 
the absence of any staggered field, for small fields is proportional to the square of the induced 
static staggered magnetic moment and the strength of magnon repulsion A. In application 
to i? 2 BaNi0 5 systems, 5(H W ) is thus expected to scale proportionately to |W Nl )] 2 , as long as 
the ordered moment is small. This linear relation was confirmed for all previously studies 
(Nd x Y 1 _ a .) 2 BaNi0 5 systemsffi 

The data obtained in the present work, particularly the more accurately measured tem- 
perature dependence of the Ni-sublattice magnetization, can be used to analyze the relation 
between A and in greater detail. In order to do so, we first have to compensate for 

the intrinsic weak T-dependence of the Haldane gap energy. This is done by replacing Ao 
by Aq(T), the temperature-dependent gap in the absence of staggered field. For A (T) we 
can take the Haldane gap measured as a function of temperature in Y2BaNiC>2, that, for 
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practical purposes, we shall approximate by the following empirical formula:i'il 

A (T) = A (0) + v /^Texp(-Ao(0)/ K T) (12) 



The coefficients A (0) = 9.21 meV and (3 = 1.0 meV 2 /K were obtained by fitting this form 
to the experimental data for Y2BaNiO"20^ This formula is used to plot 5(H W ) as a function 
of [W Nl )] 2 for all four (Nd^Yi-rr^BaNiOs systems studied (Fig. |14"D . A perfect data collapse 
is obtained. Also apparent is the slow upward curvature of this function that could not be 
clearly seen when only powder data were available.!^ 



In Ref. [L8] no quantitative predictions regarding the slope of the linear relation between 
5(H n ) and [m^ Nl ^] 2 , i.e., the coefficient A in Eq. |TT], were made. Instead, A was treated as 
an adjustable parameter. In a more refined approach A is to be determined independently 
from the properties of the Renormalization Group fixed point. This coefficient was recently 
calculated numerically^! The non-linearity in the relation between 8(H n ) and |W Nl )] 2 can 
be accounted by including higher-order terms into the 0-expansion of the Lagrangian. The 
corresponding coefficients, up to the 8th power in phi, were also recently computed.^ These 
numerical results were then utilized to calculate the staggered- moment dependence of the 
gap energy for both transverse-polarized and longitudinal Haldane excitations:0 

A \~ A ° = 0.395(m< Ni )) 2 + 0.156(m( Ni )) 4 + 0.049(m( Ni )) 6 (13) 
Ao 

A 2 - A 2 

11 ° = 1.185(m( Ni )) 2 + 0.78(m( Ni )) 4 + 0.343(m( Ni )) 6 (14) 
Ao 

These equations have no adjustable parameters and for Nd2BaNiOs systems yield the curves 



shown in solid lines in Fig. [14]. We see that our data almost perfectly agrees with the 
prediction for the transverse spin gap. The deviations at large m^ Nl ^ are expected, as the 
series in Eqs. (13|,14) are terminated at the 3rd term. Assuming that what we observe in 
our experiments are transverse gap excitations, the HC/SF model thus gives quantitatively 
correct predictions for the temperature dependence of the gap energies. 
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V. CONCLUSION 



The results of three independent series of experiments, powder neutron diffraction, and 
thermal- and cold- neutron inelastic scattering, are all consistent with our HC/SF model for 
_R 2 BaNi05 compounds at the quantitative level. The interactions between the Ni and rare 
earth subsystems in (Nd a ,Y 1 _ 2 .) 2 BaNi05 can be to a very good approximation treated at the 
MF level. This description appears to apply very well not only to static magnetic properties, 
but also to the spin dynamics. The dispersionless single-ion character of low-energy spin 
waves is a manifestation of an effective separation of the Ni-chain spin dynamics from that 
of the Nd-subsystem. 
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FIGURES 

FIG. 1. Typical measured (symbols) and calculated (solid line) powder diffraction profiles for 
in (Nd x Yi_ :E )2BaNi05 samples with x = 1 (a), x = 0.75 (b), x = 0.5 (c), and x = 0.25 (d), and the 
corresponding residuals (obs - calc). For x = 1, the full diffraction pattern (nuclear and magnetic) 
is shown. For the three other compounds, the non-magnetic scattering profile was measured just 
above the Neel temperatures and subtracted from shown scans. The ticks indicate positions of 
Bragg reflections, both nuclear (N) (in the x = 1 compound) and magnetic (M). In the magnetically 
ordered phase of the x = 0.25 compound (d), additional intensity is also observed at nuclear Bragg 
peak positions (marked with an *). 

FIG. 2. A schematic view of the magnetic structure of Nd2BaNiOs. A single crystallographic 
unit cell is shown. The magnetic moments of both Ni 2+ and Nd 3+ are confined to the (a, c) 
crystallographic plane. The magnetic propagation vector is 

FIG. 3. Measured ordered moment on the Ni (solid circles) and Nd (open circles) sites in 
(Nd^Y^^BaNiOs samples for x = 1, x = 0.75, x = 0.5, and x = 0.25, plotted against tempera- 
ture (symbols). Solid lines are calculated from the mean field model described in the text. 

FIG. 4. Measured angle formed by the Ni (solid circles) and Nd (open circles) ordered moments 
with the c crystallographic axis in (Nd^Yi-^BaNiOs samples for x = 1, x = 0.75, x = 0.5, and 
x = 0.25, plotted against temperature. Solid lines are guides to the eye. 

FIG. 5. Sample inelastic scans measured in single-crystal Nd2BaNi05 (sample A) at the 1-D 
antiferromagnetic zone-centers (open symbols) and background scans taken slightly off these posi- 
tions (solid symbols) at T = 50 K. The shaded areas indicates the scan range contaminated by a 
2ki — 3k f spurion. 

FIG. 6. Sample inelastic scans (background subtracted) measured in single-crystal Nd2BaNiOs 
(sample A) at the 1-D antiferromagnetic zone-centers (open symbols). The solid lines are fits with 
a model cross section as described in the text. 
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FIG. 7. Sample inelastic scans (background subtracted) measured in single-crystal Nd2BaNiOs 
(sample B) at the 1-D antiferromagnetic zone-centers (open symbols) at several temperatures. The 
solid lines are fits with a model cross section as described in the text. 

FIG. 8. Measured temperature dependencies of the gap energy (top) and energy-integrated 
intensity (bottom) of Ni-chain Haldane excitations in single-crystal Nd2BaNiOs. Triangles and 
circles show data obtained in the two samples studied. The solid lines are guides to the eye. 

FIG. 9. Sample inelastic scans measured in Nd2BaNiOs (sample B) at the magnetic Bragg 
position (0.5,0,1.5) using thermal (top) and cold (bottom) neutrons at T = 10 K. Apart from 
the resolution-limited feature at 4 meV energy transfer, no inelastic scattering is observed above 
background level up to 8 meV energy transfer. Solid lines are guides to the eye. 

FIG. 10. Sample inelastic scans measured in Nd2BaNiOs (sample B) at the magnetic Bragg 
zone-center (a) and zone-boundaries (b,c) illustrate the absence of any dispersion in the 4 meV 
excitation at T = 10 K. Solid lines are Gaussian fits. 

FIG. 11. Temperature evolution of the 4 mev excitation measured in Nd2BaNiOs (sample B) 
at the magnetic Bragg position (0.5,0,1.5). The solid lines represent empirical fits as described in 
the text. At T > Tjv = 48 K the excitation is totally unobservable. 

FIG. 12. Measured temperature dependence of the energy squared (open circles) and full energy 
width at half maximum (solid circles) of the 4 meV excitation in Nd2BaNiOs (sample B). The solid 
line is a theoretical curve described in the text. 

FIG. 13. Measured ordered Ni 2+ -moment in (NdzYi-^BaNiOs compounds plotted against 
the average ordered moment on the i?-sublattice (symbols, bottom axis). This dependence is 
interpreted as the staggered magnetization curve for an isolated Haldane spin chain. The estimated 
staggered exchange field acting on the Ni-subsystem is shown on the top axis. The solid line is an 
empirical fit to the experimental points, as described in the text. 
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FIG. 14. Measured increase in the Haldane gap energy in (NdzYi-^BaNiOs compounds 
relative to that in Y2BaNiOs plotted as a function of the square of the ordered staggered moment 
on the Ni-chains (symbols). For Nd2BaNiOs single-crystal data obtained in this work is shown. 



For the three other materials the data are taken from Ref. |15|. The solid lines are theoretical 
predictions for the longitudinal and transverse gaps (Ref. 
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